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Abstract 

A vast literature on the theory and phenomenology of Two-Higgs-Doublet models (2HDM) exists since 
long. However, the present situation demands a revisit of some 2HDM properties. Now that a 125 GeV 
scalar resonance has been discovered at the LHC, with its couplings to other particles showing increasing 
affinity to the Standard Model Higgs-like behavior, the 2HDM parameter space is more squeezed than ever. 
We briefly review the different parametrizations of the 2HDM potential and discuss the constraints on the 
parameter space arising from the unitarity and stability of the potential together with constraints from the 
oblique electroweak T-parameter. We also differentiate the consequences of imposing a global continuous 
U(l) symmetry on the potential from a discrete Z 2 symmetry. 


1 Introduction 


The discovery of a new boson in July 2012 by the ATLAS [1] and CMS Collaborations [ 2 ] of the CERN Large 
Hadron Collider (LHC) is undoubtedly the greatest achievement of this decade in the field of Particle Physics. 
This is most likely ‘the’ Higgs boson [3-7], the so far eluding final missing piece of the Standard Model (SM). But 
the SM has certain inadequacies. For example, it cannot account for observations like neutrino oscillations and 
dark matter. It cannot also provide adequate matter-antimatter asymmetry of the universe. These constitute 
the primary motivation to look for avenues beyond the SM, which we often call BSM scenarios. The SM relies 
on the minimal choice of a single SU( 2 ) scalar doublet acquiring a vacuum expectation value (vev) for giving 
masses to all the particles (except the neutrino) contained in the SM. One natural direction towards constructing 
BSM scenarios is to extend the SM scalar sector. In doing so, one may run into the risk of altering the tree 
level value of the precisely measured oblique electroweak parameter p (or, equivalently T). If we construct an 
SU( 2 ) x U(l) gauge theory with N number of scalar multiplets, then the general expression for the tree level 
p-parameter is [ 8 ] 
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(i) 


where and Yi denote the weak isospin and hypercharge of the i-th scalar multiplet respectively, and Vi is 
the vev acquired by the neutral component of the that multiplet. It is easy to verify that if the scalar sector 
contains only SU( 2 ) singlets (2J = 0) and doublets (2J = 1 / 2 ) with Yi — 0 and ±1 respectively, then p tree = 1 
is automatically satisfied without requiring any fine tuning among the vevs. This conforms to the experimental 
value of p, which is very close to unity [9]. In this article we restrict our discussions to the doublet extensions 
only. The simplest extension of this type is two Higgs-doublet model (2HDM) [10], which has received a lot of 
attention mainly because minimal supersymmetry relies on it. 2HDM scenarios have also been investigated to 
look for additional sources of CP violation for generating baryon asymmetry of the universe of sufficient size [ 11 ]. 
In a general 2HDM, both the scalar doublets, which we call <f>i and <f> 2 5 can couple to fermions of both types 
with X 3 = 1/2 (up-type) and — 1/2 (down-type). The up- and down-type Yukawa matrices are not, in general, 
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simultaneously diagonalizable. This will introduce flavor changing neutral currents (FCNC) mediated by the 
neutral scalars at tree level. It was shown by Glashow and Weinberg [12], and independently by Paschos [13], 
that such tree level FCNC can be avoided if fermions of a particular electric charge receive their masses from 
a single scalar doublet. This prescription can be realized by introducing a discrete or a continuous symmetry 
that apply on the scalars <Fi and <£>2 as well as on the fermions. Under the discrete symmetry one of the 
scalars is even (<£>2 —>* $ 2 ) and the other is odd (<£q —$ 1 ). There are four different possibilities for assigning 

Z 2 parities to the fermions which can avoid tree level FCNCs, i.e. ensure natural flavor conservation, so that 
Glashow-Weinberg-Pashcos theorem holds true. These correspond to the following four types of 2HDMs: 

(i) Type I: All quarks and leptons couple to only one scalar doublet <£> 2 - 

(ii) Type II: <£>2 couples to up-type quarks, while <£>i couples to down-type quarks and charged leptons (minimal 
supersymmetry conforms to this category). 

(iii) Type X (or ‘lepton specific’): 4>2 couples to all quarks, while <£>i couples to all leptons. 

(iv) Type Y (or ‘flipped’): 4>2 couples to up-type quarks and leptons, while 4>i couples to down-type quarks. 

There is also the option for preventing tree level FCNC by assuming the up- and down-type Yukawa matrices 
to be proportional to each other [14]. However, the radiative stability of the absence of FCNC couplings in 
these models is not guaranteed [15]. The Branco-Grimus-Lavoura (BGL) model [16] does on the other hand 
admit tree level FCNC couplings. But those couplings are related to the off-diagonal entries of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix and are naturally suppressed. The phenomenology of the BGL scenario 
has been studied in detail in Refs. [17,18]. A conceptually similar idea for suppressing FCNC with discrete 
symmetries was pursued in [19]. In this article, we will not elaborate any further on the FCNC issues, as we 
will not discuss the Yukawa sector of 2HDM. 

In the subsequent sections, we analyze the scalar potential, identify the physical scalar eigenstates, and reach 
the limit in which one physical scalar resembles the 125 GeV Higgs boson. For simplicity, we assume all the 
parameters in the potential to be real so that CP is manifestly conserved in the scalar sector. For discussions 
on CP violation in 2HDM scalar sector we refer the reader to Refs. [20-22], where conditions for CP viola¬ 
tion/conservation have been diagnosed in detail. With the above assumption, we derive the relations among 
the parameters of the potential that need to be satisfied to ensure that the potential is stable, i.e. it is bounded 
from below, at the weak scale. We then derive the constraints arising from the requirement of unitarity by 
studying the scattering amplitudes of 2 —>> 2 states involving the scalars and the gauge bosons. After that we 
combine the stability and unitarity constraints to impose numerical constraints on the physical scalar masses 
and other parameters. Now that one scalar has been observed around 125 GeV, constraints on the remaining 
parameter space have become more stringent. We conclude by highlighting some salient features that arise from 
the above considerations. 


2 The scalar potential 


There are two equivalent notations that are used in the literature to write the 2HDM scalar potential with a 
softly broken Z 2 symmetry (<£>1 —>• 4>i, <£>2 —>• — ^> 2 ) : 
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The bilinear terms proportional to m \ 2 in Eq. (2) or A 5 in Eq. (3) breaks the Z 2 symmetry softly. The presence of 
these soft breaking terms has implications in ensuring decoupling behavior of these models (briefly discussed in 
the concluding section). Note that when we minimize the potential of Eq. ( 2 ), the two minimization conditions 
can be used to trade and m 22 for v\ and V 2 and the potential can be cast in the form of Eq. (3). The 
connections between the parameters of Eq. (2) and Eq. (3) are given below : 


TO u = —(Aif? + X3V 2 ) ; TO22 = -(A2U2 + X3V 2 ) ; m\ 2 = y viv 2 ; ft = 2 (Ai + A 3 ) ; 

ft = 2(A2 + A 3 ) ; ft = 2 A 3 + A 4 ; ft = - —— — A 4 ; ft = ——— • (4) 


In Eq. (4), v = \/vf v 2 = 246 GeV, where v\ and V 2 are the vevs of the two doublets 4>i and 4>2 respectively. 
For most part of this article, we choose to work with the notation of Eq. (3). 

Before we proceed further, we comment on the relative status of the two parametrizations of the potential. The 
structure in parametrization 1 is more general than that in parametrization 2. It is possible to go from ‘E to 
‘2’ but not the other way. In the second one it has been assumed that both scalars receive vevs, while for the 
first this need not be the case. So the inert doublet scenario can be realized only in the first parametrization, 
e.g. in a simple illustrative scenario, when the dimensionless couplings = P 3 = ^4 = /3s = 0 , the mass mixing 
parameter m \ 2 = 0, and m 22 > 0. Then the second Higgs doublet 4>2 does not acquire any vev, and the SM 
scalar potential is recovered with the relation v 2 = v\ = —rn\ x /^\. On the other hand, If both doublets do 
indeed get vevs, there is a correspondence between the parameters in the two cases, which we have explicitly 
written down in Eq. (4). 


2.1 Physical eigenstates 

We express the scalar doublets as 


= 2 _ ( ^ w t . 

y/2 \(hi + v i) + iZi 


( 5 ) 


Then we construct the mass matrices using Eq. (3). Since we have assumed all the potential parameters to be 
real, there will be no bilinear mixing term of the form hiZj. As a result, the neutral mass eigenstates will also 
be the eigenstates of CP. For the charged sector we get the following mass matrix : 


Rma!s rSed = K 



with = — 
c 2 


u 2 

-V 1 V 2 



( 6 ) 


Diagonalizing we obtain a physical charged Higgs pair (H±) and a pair of charged Goldstones (cj^) as 
follows : 


(_ f cos/3 sin/3\ 

sin/3 cos/3y 5 


where tan/3 = V 2 /V 1 . The mass of the charged Higgs pair ( H ± ) is found to be 


2 ^4 2 

m H+ = ~2 V ■ 


(8) 
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Similarly for the pseudoscalar part one can easily find 


v: 


CP odd 


= (z x Z 2 ) -M 2 p 
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( 9 ) 


The diagonalization in CP odd sector is similar to that in the charged sector. Here we will get a physical 
pseudoscalar (A) and a neutral Goldstone (() as follows : 

C\ _ f cos p sin/A (z\ 

A I l—sin/3 cos f3 I \z 2 


( 10 ) 


The mass of the pseudoscalar is given by 


2 A 6 2 

m A = y v 1 


For the CP-even scalar part we find 


v' 


CP even 


= (hi h 2 ) 


^ 2 ) 


with, Mg = 


A s Bs 
B s C s 


where, 


As — 2 (Ai + A 3 )v\ + ~^ v 2 : 

Bs = 2(A 3 + Y)vxv 2 , 

Cs = 2(A2 + \z)v 2 + p v i ■ 

The masses of the physical eigenstates, H (heavier) and h (lighter), can be readily obtained as 


(ii) 

( 12 a) 
( 12 b) 
( 12 c) 
( 12 d) 

(13a) 
(13b) 

(14) 

As ~ C s Ai v\ - \ 2 v\ + (A,j + Ai) (>7 - v\) 

Note that there were eight parameters to start with: t?i, v 2 and 6 lambdas. We trade v\ and v 2 for v and tan/3. 
All the lambdas except A 5 can be traded for 4 physical scalar masses and a. The relations between these two 
equivalent sets of parameters are given below : 
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mH 2 
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(As + Cs) - \J(As — Cs) 2 + B 2 S 


The physical scalars are obtained by rotating the original basis by an angle a : 
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This rotation angle is defined through the following relation 
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\m\ cos 2 a + m 2 H sin 2 a — sin a cos a tan /3 (m 2 H — m\) ] — p (cot 2 /3 — l) 


(16a) 

(16b) 

(16c) 

(16d) 

(16e) 


Among these, v is already known (246 GeV) and if we assume that the lightest CP-even Higgs is what has been 
observed at the LHC, then mh is also known (125 GeV). The rest of the parameters need to be constrained 
from theoretical as well as experimental considerations. 
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Figure 1: The red, orange and yellow regions represent the 68%, 95% and 99% C.L. allowed regions, respectively, 
coming from the Higgs signal strength measurements at the LHC. The left (right) panel shows the situation for 
Type I (II) model. The star-marked points correspond to the best fit values (figures adapted from [23]). 


2.2 The alignment limit 

The alignment limit corresponds to recovering a CP-even scalar mass eigenstate with exactly the same gauge, 
Yukawa and self couplings couplings at tree level as those of the SM Higgs bosons. To start with, it is instructive 
to look at the trilinear gauge-Higgs couplings which stem from the Higgs kinetic terms: 

= l^$l| 2 + \D^2\ 2 3 y W+W^^fH + V 2 h 2 ) . (17) 

Clearly, the combination 




H° = -(v 1 h 1 +V 2 h 2 ) (18) 

V 

will have gauge couplings exactly as the SM Higgs boson and its orthogonal combination (R) will not have any 
RZZ or RWW trilinear couplings. H° also mimics the SM Higgs in Yukawa couplings. The states H° and R 
can be obtained by applying the same rotation as in the charged and pseudoscalar sectors: 

f H °\ _ ( cos ^ sin/?\ (hi\ , , 

\R) \-sinp cos P) \h 2 ) ' [ ’ 

Note that this SM-like state H° is not guaranteed to be a mass eigenstate in general. The alignment limit 

specifically implies the condition under which H° coincides with one of the CP-even physical eigenstates. To 

go to the physical basis ( H , h) from (iJ°, R ) one needs the following rotation: 

H = cos(/3 — a)H° — sin(/3 — a)R , (20a) 

h = sin(/3 — a)H° + cos(/? — a)R. (20b) 

Clearly, if we want the lightest CP-even scalar h to posses SM-like couplings, we must set sin (/3 — a) = 1, which 
is the definition of the alignment limit. Thus by going to this limit one more parameter is reduced. 

Now we come to the important question of how crucial the alignment limit is in the context of current LHC 
Higgs data. Many global fit results in view of the recent data can be found in the literature [23-28]. In Fig. 1 
we display the results of a recent analysis [23]. The orange part represents the 95% C.L. allowed region from 
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Figure 2 : Current measurements of the Higgs signal strengths into different channels by the ATLAS [30] and 
CMS [31] Collaborations. 


measurements of the Higgs signal strengths in various channels (for the latter, see Fig. 2 ). Since the data is 
compatible with the SM predictions, the alignment limit is preferred. The horizontal widths of the allowed 
regions reflect the present accuracy of measurements. In Ref. [28], it has been shown how this width will shrink 
if future measurements continue to agree more with the SM predictions with greater accuracy, thus pushing us 
closer to the alignment limit. The SM alignment limit has also been motivated by employing different global 
symmetries of the scalar potential [29]. 


3 Stability constraints 


To discuss the stability of the potential it is convenient to work with the notation of Eq. (2). Here we derive 
the constraints on parameters Pi such that the scalar potential V is bounded from below in any direction in the 
field space [32,33]. It is sufficient to examine the quartic terms of the scalar potential (which we denote by Vf) 
because only this part of the potential will be dominant for large values of the field components of Tq and <I> 2 - 
We define a = Tjdq. 6 = <f> 1&2, c = Re d = Im <£>{<!>2 and note that 

ab > c 2 + d 2 . (21) 

Using these definitions we can rewrite the quartic part of the scalar potential as follows [34]: 

V 4 = - (^/p[a — \f~P 2 b'j + ^3 + \/P 1 P 2 ^ (ab — c 2 — d 2 ) + 2 ^3 + @4 + y/Pip 2 ^j c 2 

+ (Re /? 5 - /? 3 - ^ - y/?i/? 2 ) (c 2 - d 2 ) - 2 cd Im /? 5 . (22) 

Although we assume all the potential parameters to be real for our phenomenological studies, our arguments 
on stability do not depend on whether ^5 is real or complex. We must ensure that V 4 never becomes infinitely 
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negative in any direction of the field space, i.e ., for any choice of 8 independent field parameters (4 of 4>i and 4 
of 4> 2 ). Since 4>i and 4 > 2 are two component column matrices, it is possible to choose arbitrary nonzero values 
for a and b even when we make c = d = 0. But if a and/or b becomes zero, then c = d = 0 automatically. 
Keeping these in mind, we now proceed to find the stability constraints, i.e. the conditions under which the 
potential is bounded from below [32]. 


(i) Consider the field direction b 

= 0 (and 

therefore c 

= d = 

= 0 ) and a — 

oc; then V4 = 

= (ft /2 )a 2 . 

So, V4 is 

not largely negative requires 










0 i 

> 0 . 




(23) 

(ii) Consider the field direction a 

= 0 (and 

therefore c 

— d - 

= 0 ) and b -+ 

oc; then V4 = 

= (Mb 2 . 

So, V4 is 

not largely negative requires 










02 

> 0 . 




(24) 


(Hi) Consider the field direction along which a = yJfc/Pib (so that the first term in Eq. ( 22 ) vanishes) and 
c = d = 0. In addition to this we go to large field values in that direction, i.e., a, b -+ oo. Then, 
V 4 = (03 + \[0\ 02 )ab. Now, as a, b > 0 by definition, the condition for the potential not to hit (— 00 ) 
becomes 

ft + vTft > 0. (25) 


(iv) 


Again consider the field direction in which a = \ff+W\b along with ab = c 2 + d 2 . Along this direction, 
V4 is of the form 


V 4 = Pc 2 + 2 Qcd + Rd 2 , (26a) 

where, P = Re 0$ + A, (26b) 

Q = -Im 0 s , (26c) 

R = -Re 05 + A, (26d) 

with, S = 03 + 04 + \j0\0 2 • (26e) 

Since c and d are still arbitrary, by choosing d = 0 , c -+ 00 and c = 0 , d -+ 00 successively, we require 

P = Re 05 + S > 0, (27a) 

R = -Re 05 +S> 0 , (27b) 

and hence, S > 0 . (27c) 


To have another condition, let us recast Eq. (26a) into the following form : 

V 4 = P(c+|d) +(R-^Pjd 2 . (28) 

We can now choose a direction along which c— —Q/P d with d -+ oc so that we have the following condition : 

O 2 

R-^> 0 => PR>Q 2 . (29) 

For the last step, remember that P > 0 (Eq. (27a)) so that we can multiply both sides by P without flipping 
the inequality sign. After substituting for P, Q and R we get from Eq. (29) : 

S 2 - (Re k ) 2 > (Im kf ^S 2 > |/3 5 | 2 , (30a) 

S > |As|, (30b) 

where, in the last step we have used the fact that S > 0 (Eq. (27c)). Since \0$\ > +05 , 0, Eq. (30b) puts a 
stronger constrain on S than Eq. (27). Therefore, substituting for S, Eq. (30b) becomes 

03 T 04 T \J 0102 > \05 | • (31) 


7 








We now collect Eqs. (23), (24), (25) and (31) together and, using Eq. (4), express them in terms of lambdas for 
later use : 


Ai + A 3 > 0 , 

(32a) 

^2 + A 3 > 0 , 

(32b) 

( 2 A 3 + A 4 ) + 2a/(Ai + A 3 /A 2 + A 3 ) > 0 , 

(32c) 

2 As+ A 5 + A 8 l As 2 A «l + 2,/fA, + A 3 )(A 2 + A 3 ) > 0. 

(32d) 


The above conditions are both necessary and sufficient for ensuring the stability of the electroweak vacuum. 
This has also been shown through rigorous analysis in [35,36]. In addition to these stability conditions, one 
also needs to ensure the positivity of the physical scalar masses [33,37]. Additionally, one might also wish to 
ensure that the minima is indeed the global minimum. The condition for the latter is given by [38] 



4 Unitarity constraints 


In this section we study the energy growth of scattering amplitudes involving the scalar states. Any scattering 
amplitude can be expanded in terms of the partial waves as follows : 

oo 

M{6) = 167T a^(2£ + l)P/cos 0 ), (34) 

£=0 

where, 6 is the scattering angle and Pi(x) is the Legendre polynomial of order £. The prescription is as follows: 
once we calculate the Feynman amplitude of a certain 2—^2 scattering process, each of the partial wave 
amplitude (a/, in Eq. (34), can be extracted by using the orthonormality of the Legendre polynomials. In 
the context of SM, the pioneering work has been done by Lee, Quigg and Thacker (LQT) [39]. They have 
analyzed several two body scatterings involving longitudinal gauge bosons and physical Higgs in the SM. All 
such scattering amplitudes are proportional to Higgs quartic coupling in the high energy limit. The £ = 0 partial 
wave amplitude (ao) is then extracted from these amplitudes and cast in the form of what is called an S-matrix 
having different two-body states as rows and columns. The largest eigenvalue of this matrix is bounded by the 
unitarity constraint, |ao| < 1. This restricts the quartic Higgs self coupling and therefore the Higgs mass to a 
maximum value. 

The procedure has been extended to the case of a 2HDM scalar potential [40-44]. Here also same types of two 
body scattering channels are considered. Thanks to the equivalence theorem [45,46], we can use unphysical 
Higgses instead of actual longitudinal components of the gauge bosons when considering the high energy limit. 
The diagrams containing trilinear vertices will be suppressed by a factor of E 2 coming from the intermediate 
propagator. Thus they do not contribute at high energies, and only the quartic couplings contribute. Since we 
are interested only in the eigenvalues of the S-matrix, we may for convenience with the original fields of Eq. (3) 
instead of the physical mass eigenstates. 

To provide clarity, let us outline the method of obtaining the constraints. As already argued, only the dimen¬ 
sionless quartic couplings will contribute to the amplitudes at high energies. As a result, only £ = 0 partial 
amplitude (ao) will receive nonzero contribution from the leading order terms in the scattering amplitudes. The 
task is to find the expressions of ao for every possible 2 —)> 2 scattering process and cast them in the form of an 
S-matrix which is constructed by taking the different two-body channels as rows and columns. Unitarity will 
restrict the magnitude of each of the eigenvalues of this S-matrix to lie below unity. 

First we identify all the possible two-particle channels. These two-particle states are made of the fields hk 

and Zk corresponding to the parametrization of Eq. (5). For our calculation, we consider neutral combina¬ 
tions out of two-particle states (e.g., wf wj , hihj, ZiZj , hiZj) and singly charged two-particle states (e.g., 
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hj , Zj ). In general, if we have n- number of doublets fa (k = 1,..., n) in nHDM scenarios, there will be 

(3 n 2 +n)-number of neutral and 2 n 2 -mimber of charged two-particle states. Clearly, the dimensions of S-matries 
formed out of these two-particle states will be a (3n 2 + n) x (3n 2 + n) and 2 n 2 x 2 n 2 for the neutral and charged 
cases respectively. The eigenvalues of these matrices should be bounded by the unitarity constraint. 


The neutral channel S-matrix for 2HDM is a 14 x 14 matrix with the following two-particle states as rows and 
columns : 


Wi Wi 


tin 


WA W. 


2 5 


1 /1 


Wi W. 


2 5 


on 


WA W 


h\hi zizi h 2 h 2 z 2 z 2 

75 ’ vr 75 5 


hiz 2 , h 2 z\, z\z 2 , h\h 2 , h\Z\, h 2 z 2 . 


(35) 


The factor of 1/75 associated with the identical particle states arises due to Bose symmetry. In the most 
general case, finding the eigenvalues of the 14 x 14 matrix would be a tedious job. But the potential of Eq. (3) 
contains some obvious symmetries in its quartic terms. These symmetries will allow us to decompose the full 
matrix in smaller blocks. Now, each term in the quartic part of the potential always contains even number of 
indices (1 or 2 ). Consequently a state x\y\ or x 2 y 2 will always scatter into x\y\ or x 2 y 2 but not into x±y 2 or 
x 2 yi and vice versa. Furthermore, CP symmetry is conserved. This implies that a neutral combination hihj or 
ZiZj will never go into hiZj. Keeping these facts in mind we can now decompose the S-matrix in the neutral 
sector into smaller blocks as follows : 


7*^ TV 7x6 

M. n = [ 0 

0 


The submatrices are given below : 


(A4jv )6X6 


* 1*1 

V2 

h^hi 

* 2*2 
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fa -2 h 2 
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\/2 (As + 4f) 
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A 3 + 4f 
a 3 + 4& 


h-^ 

v / 2(Ai + A 3 ) 
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(Ai + A 3 ) 
3(Ai + A 3 ) 
A 3 + 4f 
A 3 + ^ 


(36) 
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V* 
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75 (a 2 + a 3 ) 
a 3 + if 

a 3 + if 
3(A 2 + A 3 ) 

(a 2 + A 3 ) 
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75 (a 3 + %)\ 
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a 3 + if 

A 3 + 

(a 2 + A 3 ) 
3(A 2 + A 3 ) 


/37a) 


h 1 z 1 


^ 2*2 


(M 


N J 6 X 6 — 




(A#) 

h i z i ( 2(Ai 

2X2 - 1 a 5 

h 2*2 \ 

+ A 3 ) Aa 

T* 2(A 2 

2 Ae 4 
+ As)b 





'i 

v 2 

W 2 W 1 

^1*2 

^2*1 

* 1*2 

h 

1^2 

w l w 2 

( 2A 3 + 

A 5 +^6 

2 

As — A6 

i(A 4 -A 6 ) 

i(A 4 -A 6 ) 

A 5 —a 4 

2 

A s 

^ \ 

w i w i 

A 5 - 

- Xq 

2A 3 + As + Ae i 

j(A 4 — Ae) 

-§(A 4 -A 6 ) 

a 5 —a 4 

2 

A 5 

-A 4 

2 

h\z 2 

|(A4- 

- Ae) 

-|(a 4 -a 6 ) 

2A 3 + Ae 

A 5—^6 

2 

0 


0 

h 2 z 1 

-i(A 4 

— Ae) 

i(A 4 - Ab) 

A 5 —A 6 

2 

2A 3 + Ae 

0 


0 

*1*2 

^ 5 - 

2 

a 4 

a 5 —a 4 

2 

0 

0 

2A 3 + A 5 

A 5 

-A 6 

2 

h 1 /i 2 

L A 5 - 

\ 2 

a 4 

A 5 —a 4 

2 

0 

0 

A 5 — A 6 

2 

2A 3 

+ A 5 / 


(37b) 


(37c) 


The same exercise can be repeated for the charged two-particle state combinations. With the singly charged 
state combinations, it will be a 8 x 8 matrix which will take the following block diagonal form : 


\a _ ( 7 x 4 0 

MC - ( 0 (Mg) 4 y4 


(38) 
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The submatrices are given below : 


(-A^c )4x4 


(A / fJ?)4x4 


h2W^ 

Z\Wf 

h2W + 

Z!W + 
Z 2 w^ 


h\w+ 

h 2 W^ 

Z\W^ 

Z2W2 

/ 2 (Ai + A3) 

X5 — X4 

2 

0 

-|(A 4 - A 6 )\ 

A5 — A4 

2 

2(A2 + A3) 

-|(a 4 -a 6 ) 

0 

0 

|(a 4 -a 6 ) 

2 (Ai + A3) 

A5 —A4 

2 

V |(A 4 — Ae) 

0 

A5—A4 

2 

2 (A 2 ~b A3) / 

h\W 2 

h 2 W + 

Z1W2 

z 2 wf 

( 2A3 + A4 

X5 — X4 

2 

0 

5 (A4 - A 6 ) \ 

A5 —A4 

2 

2A3 + A4 

|(a 4 -a 6 ) 

0 

0 

— §(A 4 -A 6 ) 

2A3 + A4 

A5 —A4 

2 

\-|(a 4 -a 6 ) 

0 

A5 —A4 

2 

2A3 + A4 / 


The eigenvalues for these matrices are given by 


• (-M“) 6 x 6 : af, af, af. 

• (Mk) 2 x 2 ■ af. 

• {M^)q x q • &i, ^ 2 , ^47 ^ 5 , with 65 twofold degenerate. 

• (-Mc)4x 4 : a±, «3 • 

• (VtJ?)4x4 : b‘2, b 4 , 65 , & 6 - 


We also display the explicit expressions for these eigenvalues : 

CL ] 1 = 3 (Ai + A2 + 2A3) d= 

= (Ai + A2 + 2A3) zb ]J (Ai — A2) 2 + - (2A4 — A5 — Ae) , 

a t = (Ai + A2 + 2A3) zb ]J (Ai — A2) 2 + - (A5 — Ae) , 

1 5 

bi = 2A3 — A4 — -A5 + —Xq , 

62 = 2A3 + A4 — -A5 + —Xq , 

5 1 

63 = 2A3 — A4 + -A5 — -Xq , 

64 = 2A3 + A4 + -A5 — -Xq , 

^5 = 2A3 + -A5 + -Xq , 

b§ = 2(A3 + A4) — -A5 — -A6 • 

Each of the above eigenvalues will be bounded from the unitarity constraint as 

\af I, \bi\ < 167r. 


/ 9 (Ai — A2) 2 + 4A3 + A4 + 


A; 


+ A 6 \ " 


(39a) 


(39b) 


(40a) 

(40b) 

(40c) 

(40d) 

(40e) 

(40f) 

(40g) 

(40h) 

(40i) 

(41) 


10 




Log|^iiinp> 



0 300 600 900 1200 1500 

m A *CeV> 


Figure 3: Allowed region (shown by scattered points) from unitarity and stability for exact Z 2 symmetry (\$ = 0/. 


5 Numerical constraints on the scalar masses 


We now investigate the implications of the above conditions on the physical scalar masses, especially the non¬ 
standard ones. Fig. 3 shows the region allowed by the combined constraints from unitarity and the boundedness 
of the potential for A 5 = 0, i.e ., exact Z 2 symmetry. Two noteworthy features emerge [47] : 

( 1 ) From the left panel of Fig. 3, one can read the upper and lower limits on tan/3 as 1/8 < tan/3 < 8 . 

(ii) The upper limits on nonstandard scalars masses are given by: m#, m#+ < 1 TeV. 

To understand the origin of the above limits we look into the eigenvalues of Eq. (40). The first two constraints 
for boundedness in Eq. (32) can be combined into 

Ai + A 2 + 2 A 3 > 0 . (42) 

Together with the condition \af\ < I 67 r, this implies 

0 < Ai + A 2 + 2 A 3 < —— , (43) 

o 

=> 0 < (m 2 H - ~ X^v 2 ^ (tan 2 /3 + cot 2 /3) + 2 m\ < , (44) 

where the last expression is obtained from the previous one by using Eq. (16) in the alignment limit. Since the 
heavier CP-even Higgs mass ra# > 125 GeV, a limit on tan /3 (as well as cot /3) is obtained when A 5 = 0. Since 
the minimum value of (tan 2 /3 + cot 2 /3) is 2 when tan /3 = 1 , the maximum possible value of ra# is obtained for 
tan/3 = 1. Eq. (44) thus explains the tan/3 dependent bound on ra# as depicted in the left panel of Fig. 3. 

To obtain restrictions on the individual masses, tra and m H +, we use the following inequalities : 

1 167TU 2 

\h ~ h\ = 3 |A 6 - A 5 | < 327T => \m 2 A --X 5 v 2 \<—-—, (45a) 

1 167T7 ^ 2 

be ~ h\ = 3|A 4 - A 5 | < => \m 2 H+ --X 5 v 2 \ <—-—. (45b) 

Using Eqs. (45a) and (45b) we put limits on tua and ra#+, respectively, when A 5 = 0. Additionally, due to the 

inequality 

1 Gttv 2 

\bi -b 6 1 = 3|A 6 - A 4 | < 327r => \m 2 A - m 2 H+ \ < —-— , (46) 
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Figure 4: Relaxation of the constraints on tan/3 for nonzero A 5 . The right panel should be compared with the 
left panel of Fig. 3. 
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Figure 5: Effect of nonzero A 5 on the constraints on the nonstandard masses. 


we expect the splitting between and ra#+ to be always restricted. It is also interesting to note that the 
conclusions obtained from Eqs. (45a), (45b) and (46) do not depend on the imposition of the alignment condition. 

From Eq. (44) we also observe that the allowed space for tan (3 is further squeezed if A 5 < 0, but the bound is 
relaxed if A 5 >0. This feature emerges from the left panel of Fig. 4. In addition to this, we can see from Eq. 
(44) that m 2 H must be close to I/ 2 A 5 U 2 if tan/3 moderately deviates from unity. This feature is reflected by the 
horizontal tail in the right panel of Fig. 4 on both sides of the peak. The width of the tail is a result of the 
variation of A 5 in the range [-15, 15]. 

From Eqs. (45a), (45b) and (44) we note that the upper bounds on the nonstandard scalar masses will be 
relaxed for A 5 > 0 and get tighter for A 5 < 0. Fig. 5 reflects these features. 

An interesting alternative arises if instead of Z 2 one imposes an U(l) symmetry under which <f>i —>> <f>i and 
<f >2 e m 4 > 2 . This U(l) symmetry needs to be softly broken to forbid the appearance of an exactly massless 
pseudoscalar. This U(l) symmetry, in the quartic terms, implies /3s = 0 in Eq. (2) or A 5 = A 6 in Eq. (3). The 
constraints on the scalar masses imposed by the stability and unitarity conditions in Eqs. (32) and (40) have 
been plotted in Fig. 6 for tan/3 = 1, 5 and 10 by performing random scan over all nonstandard scalar masses. 

The following salient features emerge from the plots. 

(i) There is a correlation between and ra# which gets stronger for larger values of tan/3. They become 
nearly degenerate once tan/3 > 10. To understand this, we observe that Eq. (44) for A 5 = A 6 reduces to 

‘-KO'j 

0 < {jn 2 H — 777 / 4 )(tan 2 P + cot 2 fi) + ^ m h — —^— • (47) 

Clearly, for tan /3 away from unity, H and A are nearly degenerate. 
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Figure 6: 2 HDM potential with softly broken U(l) symmetry: regions allowed in mn-rriA, mH-m H + and m A - 
m H + planes from unitarity and stability (red points), and from T-parameter (black points), for three choices of 
tan [3. The plots have been taken from [48], where m H + > 100 GeV was assumed to respect LEP direct search 
bound [49]. 
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(ii) There is a similar correlation between uih and m H +, but unlike the previous point, without any depen¬ 
dence on tan (3. This can again be seen from the inequalities of Eqs. (45a) and (45b) keeping in mind that 
now rn\ = (l/2)\$v 2 . 

(Hi) The unitarity conditions essentially apply on the difference of the nonstandard squared masses. Any 
individual mass can be arbitrarily large without affecting the unitarity conditions. This conclusion crucially 
depends on the existence of a U(l) symmetry and its soft breaking term in the potential. When the 
symmetry of the potential is only a discrete Z 2 , considerations of unitarity do restrict the individual 
nonstandard masses as has already been shown. 


(■ iv ) The splitting between the heavy scalar masses is also constrained by the oblique electroweak T-parameter, 
whose expression in the decoupling limit is given by [50,51] 

1 


T = 


167T sin 2 0 - 


■jMfr L- 


F(m 2 H+ ,m 2 H ) + F(m 2 H+ ,m 2 A ) - F(m 2 H ,m 2 A ) 


(48) 


with 


Z7V \ x + y °°y w / ^ 
F(x,y) = —-In (x/y ). 

2 x-y 

The new physics part in the T-parameter is given by [52] 

T = 0.05 ±0.12. 


(49) 

(50) 


To provide intuition into the constraints from the T-parameter, we assume uih = UC 4 , which is anyway 
dictated by the unitarity constraints for tan /3 somewhat away from unity. It then follows from Eq. (48) 
that the splitting between m H + and ra# is approximately 50GeV, for \m H + — ra#| <C It 

follows from Fig. 6 that the constraints from the T-parameter are stronger than those from unitarity and 
stability. 

On the other hand, for tan /3 = 1, unitarity and stability do not compel mn and ttia to be very close. Then 
the T-parameter cannot give any definitive constraints in the planes of the nonstandard scalar masses, 
unlike the unitarity and stability constraints. For this reason, we have shown only the unitarity/stability 
constraints in Fig. 6 for tan/3 = 1. 

(v) For moderate or large tan/3, the unitarity and stability constraints together with the T-parameter con¬ 
straints imply that all three heavy scalar states are nearly degenerate in the alignment limit. 


For completeness, we also comment on the high scale validity of 2HDMs [53]. After the Higgs discovery, it is well 
known that the SM perhaps needs to be augmented by new physics at energies beyond 10 8 -10 10 GeV [54,55]. 
The reason behind this is the following: the Higgs boson mass ~ 125 GeV turns out to be a little ‘smaller’ than 
what could have made the SM ultraviolet safe up to Planck scale. More precisely, in the tussle between the 
top-Yukawa interaction (which tries to pull down the scalar quartic coupling) and the quartic self interaction 
(which tries to push itself up) during the course of RG evolution, the top-Yukawa takes an upper hand and 
pulls the SM quartic coupling to negative values well below the Planck scale. One might expect the situation 
to improve in a 2HDM because of the presence of additional quartic couplings. But the 2HDM potential of 
Eqs. (2) and (3), without the soft breaking term, still fails to maintain stability all the way up to the Planck 
scale [56-58]. Even in the presence of a soft breaking parameter, in the alignment limit, a lower bound on tan /3 
(> 3) [57,58] is obtained from the requirement that the effect of the top-Yukawa (= y/2m t /(vsinj3) in 2HDM 
vis-a-vis y/2m t /v in SM) is sufficiently diluted to maintain the required stability. Moreover, certain correlations 
between the nonstandard masses and the soft breaking parameter need to be maintained. 


6 Conclusions and outlook 


Based on data on the scalar resonance observed at the LHC, two important conclusions have been drawn: (i) 
its mass is approximately 125 GeV and (ii) its couplings to gauge bosons and fermions are SM-like. In the 
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2HDM framework, the latter observation pushes us to the alignment limit. Put these two conditions together, 
the 2 HDM parameter space is more constrained then ever. Some important observations in this context are the 
following: 

• When the potential has a global U(l) symmetry, rather than a discrete Z 2 , as well as a soft breaking term, 
unitarity restricts the mass-squared differences of the nonstandard scalars. So the individual nonstandard 
masses, like mn^TUA or m#+, can grow very large without necessarily violating unitarity. 

• Strictly when both doublets receive vevs and the potential has an exact discrete or a global continuous 
symmetry, the model exhibits features of nondecoupling. For example, observables do not necessarily 
reduce to their SM values even when nonstandard scalars are too heavy. In [59], we have demonstrated 
this behavior in the context of h 77 , where in the limit of an exact discrete symmetry, the Higgs signal 
strength in the diphoton channel at the LHC, given by /i 77 , was shown to retain nonstandard effects even 
when the charged Higgs mass is pushed to extremely large values. Employing the soft breaking term, it is 
possible to ensure decoupling. If the symmetry to start with is a discrete Z 2 , one cannot avoid fine-tuning 
between the charged Higgs mass and the soft breaking parameter to reach the decoupling limit. On the 
other hand, if the starting symmetry is a global U(l), the soft breaking parameter gets related to the 
pseudoscalar mass. In this case, the combined constraints from unitarity and the T-parameter naturally 
lead to the decoupling limit without any need for fine-tuning. For details, we refer the readers to Ref. [59]. 

Acknowledgments: DD thanks Department of Atomic Energy, India, for financial support. 
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